The Global Kinematics of the Globular Cluster M92 



G. A. Drukier 

Department of Astronomy, Yale University, New Haven, CT 
drukierOastro . yale . edu 
^ I H. N. Cohn and P. M. Lugger 

o , 

£S) \ Department of Astronomy, Indiana University, Bloomington, IN 

> 

O . cohn, lugger@astro.indiana.edu 

S. D. Slavin 

Dept. of Chemistry and Physics, Purdue University Calumet, Hammond, IN 



> 

slavinOcalumet . purdue . edu 

■ 

R. C. Berrington 

■ Dept. of Physics and Astronomy, University of Wyoming, Laramie, WY 

o 

rberringOuwyo . edu 

Q ■ and 

B. W. Murphy 

Department of Physics and Astronomy, Butler University, Indianapolis, IN 

X" 

$— i ' bmurphy@butler.edu 

ABSTRACT 

We report the determination of high-accuracy radial velocities for 299 members of the 
globular cluster M92 using the Hydra multi-object spectrograph on the WIYN telescope. 
We have concentrated on stars outside of the central region of the cluster, located up to 
14'4 from the cluster center. Candidate members were selected for spectroscopy based 
on a photometric metallicity index determined from 3-band Washington photometry, 
also obtained with the WIYN telescope. The median error in the velocities is 0.35 
kms -1 . We find the heliocentric radial velocity of the cluster to be —121.2 ±0.3 kms . 

We have used an improved Bayesian analysis to determine the velocity dispersion 
profile of M92. The most probable profile is a cored power-law with a scale radius of 2', 
velocity dispersion at 1' of 6.3 kms -1 and outer power-law with slope —0.6. We have 
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also reanalyzed the M15 radial velocities of iDrukier et al.l (|1998l ) and find that a pure 



power-law with a 1' velocity dispersion of 8 kms" 1 and slope —0.5, and the combination 
of a power-law with slope —0.4 and scale of 7.5 kms -1 inside 9' and a dispersion of 4 
kms -1 outside, are equally likely. In both clusters there is evidence that the samples 
include escaping stars. We present results from a GRAPE-based N-body simulation 
of an isolated cluster that demonstrates this effect. We suggest additional tests to 
determine the relative importance of tidal heating and stellar ejection for establishing 
the velocity field in globular cluster halos. 

Subject headings: globular clusters: individual (M92) — globular clusters: individual 
(M15) — methods: statistical 



1. Introduction 

In globular clusters, the interplay between two-body relaxation and external tidal stresses is 



and eventual dissolution of clu s ters in the Galactic tidal field (e.g. 



Combes. Leon. &: Mevlan 



1999 



Takahashi fc Portegies Zwartll2000l : 



Gnedin. Lee. 


& Ostriker 


1999; 


Dehnen et al. 


20041). To 


inves- 



tigate these issues, we have been carrying out a program to determine the global velocity fields 
of globular clusters using the Hydra fiber-fed, multi-object spectrograph on the WIYN telescop^j]. 
With its 1° diameter field and echelle grating, Hydra is well suited for determining high-accuracy 
radial velocities of stars out t o the tidal radii of clusters. Our approach is an important complement 
to Fabry- Perot imaging (e.g. iGebhardt et al J 1 19941 . l2000l ). which can only be used efficiently in the 
central 1 — 2'. 



In a previous paper (jDrukier et al.lll998l . hereafter Paper I) we reported a new global velocity- 
dispersion profile for the prototypical collapsed-core cluster M15. We found a clear indication of a 
flattening and possible rise of the profile in the outer part of the cluster. In contrast, our anisotropic 
Fokker-Planck simul ations of isolated clusters show a smoothly declining velocity-dispersion profile 
(jDrukier et al.lll999l ). We inter preted our observation s of M15 as evidence for heating of the cluster 
halo by the Galactic tidal field. iJohnston et al.1 (|1999l ) subsequently showed that a Galactic satellite 
system undergoing tidal heating will show a break in the slope of the velocity dispersion profile at the 
radius at which unbound stars begin to dominate. However, this leaves open the question of whether 
two-body relaxation in an isolated cluster could produce a similar effect through strong scattering 
of stars from the central region into unbound orbits . Neit her our Fokker-Planck simulations nor 
the self-consistent field simulations of lJohnston et al.l (|1999l ) include strong scattering. We address 
this question in the present study through direct N-body integrations, which include all scattering 



The WIYN Observatory is a joint facility of the University of Wisconsin, Indiana University, Yale University, 
and the National Optical Astronomy Observatories. 
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mechanisms. 



We chose M92 (NGC 6341) for the second cluster in our survey, as a "normal" comparison 
to the collapsed-core cluster M15. The properties of these two clusters are compared in Table [TJ 
Unlike M15, M92 has a well-resolved core. While its cen tral luminos ity density is in the upper 
third of the distribution for all Galactic globular clusters (|Harrislll996l ). it is nevertheless at least 
an order of magnitude lower than that of M15. Both clusters are brighter than the median of the 
absolute magnitude distribution; M92 is in the upper 30%, while M15 is in the upper decile. Both 
are among the very lowest metallicity clusters. Both are located approximately 10 kpc from the 
Galactic center and nearly 5 kpc from the Galac tic plane. Both are among the oldest group of 
clusters (jChabover. Demarque. k, Saraiedinilll996l ). 

The radial velocity of M92 (-121.2 ±0.3 kms" 1 ) and its very low metallicity ([Fe/H]= -2.29) 
clearly distinguish its stars from the bulk of the field stars in the same region of the sky. Nev- 
ertheless, selecting likely cluster members for spectroscopy becomes chall enging in the outer halo 
where the vast majority of stars are nonmembers. As we pointed out in Paper II . selecting stars 
based on their position in a color-magnitude diagram (CMD) alone is not an efficient means of 
finding potential members. As we discuss in th e next section, we turned to a metallicity-sensitive 
color-color diagram in the Washington system (|Canternalll976l ) to aid our candidate selection. 



There have been a number of previous velocity studies of M92. The most extensive is that of 
Lupton. Gunn. fc Griffin! ()1985l ). who observed 49 stars in M92 usi ng the 5 m Ha l e tele scope, but, 
unfortunately, the individual velocities have not been published. iLupton et al.l (|1985l ) find little 
sign of anisotropy and only a weak {2a at best) rotational signal. The disper sion profile i t self i s 
fairly flat, but the statistics in the outer part of the cluster are particularly poor. iBeers et al.l (jl990) 
observed o nly 7 stars and the unc ertainty in the velocities is 7 kms -1 . While the present study was 
underway. ISoderberg et al.l (119991 ) reported high-precision measurements for 35 stars in M92, based 
on Hydra spectroscopy with the KPNO 4m and the WIYN telescope in the region about the Na D 
line. Thus, while some previous radial velocity information is available for M92, it is not sufficient 
for a determination of the global velocity dispersion profile, particularly in the cluster halo. Good 
quality proper mo tions for M92 stars are available for of order 700 stars located out to 19' from 
the cluster center ()Reeslll992l ; iTucholke. Scholz. fc Broschelll996l ). 



As in iPaper J . we carried out our observations on the 3.5 m WIYN telescope. We used the 



Hydra fiber-fed, multi-object spectrograph to obtain high-resolution spectra. This has resulted in 
a homogeneous set of over 300 M92 stars with new, high-accuracy, velocity measurements. 

In the absence of a full dynamical study, data sets such as these are customarily analyzed 
by binning the stars in radius and computing the dispersion in their velocities. In this way, an 
approximation to the dispersion profile is constructed. In this paper, we use the methods of 
Bayesian statistics to move beyond this in three ways. First, besides the traditional bins, we also 
examine several continuous dispersion profiles as described in §3.21 Second, our methods allow us 
to assess which of the several candidate profiles is best supported by the data. Third, by looking 
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at the posterior probability distributions for the parameters we can determine the extent to which 
they deviate from Gaussian distributions. 

We discuss our observations in the following section. It particular, we review our new pho- 
tometric selection technique using the Washington system. In £}3]we use the Bayesian techniques 
descr ibed in ^ 3] to analyze the velocity distributions of both the new M92 data and the M15 data 



from iPaper It We then do a similar analysis of an 8k N-body model to investigate the effects of 
unbound stars on the dispersion analysis. We conclude with a summary of the main results of this 
paper and a suggestion of some directions for further progress in investigating the questions raised 
by these data sets. 



2. Observations 

2.1. Photometry and Candidate Selection 

In order to observe with Hydra, it is necessary to prepare a list of target stars with accurate 
positions. For this project, the targets are those stars found on the giant branch which are likely 
to be cluster members. Our candidate list comes from three sources. For our 1996 May ru n and 



the fir st night of the 1996 June run, our candidates all came from the proper motion list of iRees 



(|1992l ). We selected stars appearing on the V vs. B — V CMD, but did not select on the basis of 
the assigned membership probabilities. For the bulk of the 1996 June run, our stars were selected 
from a list of stars found on KPNO Burrell Schmidt CCD images kindly supplied by H. Morrison. 
The images had been taken in the Washington system, which we discuss further below. Owing to 
problems with obtaining sufficiently accurate photometry from these images, we selected stars only 
on the basis of their CMD position, rather than using the color-color plot which provides metallicity 
information. We observed 444 stars in 1996 and obtained useful velocities for 427 of these. Only 
192 stars (45%) proved to be members. 



As we discussed in lPaper J . selecting stars on the basis of a CMD alone is not the most efficient 
method of finding likely members. Depending on the distance and Galactic latitude of the cluster, 
contamination by Galactic dwarfs with the colors of the cluster giant branch can be significant. 
This is a particular problem in the outer part of a globular cluster where the members are sparse 
on the sky and are outnumbered by non-members even at high Galactic latitudes. A better method 
is to take advantage of the metallicity difference between a globular cluster and the Galactic field. 



The Washington system (|Canternalll976l ) is designed to determine stellar metallicities along the 
giant branch. It consists of four broad-band filters, C, M, T\, and T<i- The last is equivalent to Kron- 
Cousins /. The T\ — T<i color was originally intended to be the temperature index, with C — M giving 
the metallicity. The C band was selected to occupy the region around 400 nm occupied by the Ca H 
and K doublet, the G band, and plentiful CN bands in these cool stars. The M band is in a spectra l 



region centered near 510 nm which is less affected by metallicity. iGeisler. Claria. &: Minnitil (|199ll ) 
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have shown that the optimal diagram for metallicity determination is the C — M vs. M—T2 color- 
color diagram. This gives us the maximum metallicity sensitivity and only requires observation in 
three broad-band colors. 

In 1997 May, we observed M92 using the imager on the WIYN telescope. We observed a 
40-field mosaic covering a 15' radius centered on the cluster, with a few fields extending out to 20'. 
We observed eac h field with th e C, M, and T2 filters. We also observed using the DDO 51 filter 
as suggested by iGeislerl ()1984l ) . This should have allowed us to separate dwarfs from giants on 
the basis of the luminosity effects in the Mg b triplet. The DDO 51 photometry lacked sufficient 
precision to permit this, possibly due error propagation in transferring the zero-point across the 
large mosaic. The DDO 51 frames were used for the astrometry described below. 

Each frame was debiased and flat-fie lded in the u s ual manner, and photo metry was done on 
each using DAOPHOT II and ALLSTAR (jstetsorJll987l : lstetson k Harrislll988h . Since the imaging 



was not done under photometric conditions, all of the photometry was put on a common zero-point 
by matching stars in overlapping regions between frames. In a subsequent, photometric, run in 
1997 October, we observe d one fie l d in M 92 as well as a region in standard field SA114 containing 
three standard stars from IGeislerl (|1996l ). This allowed us to put our C, M, and T2 observations 
onto the standar d system. Having only th e one standard image, we used the first order extinction 
coefficients from lHarris fc Canternal ()1979l ). Our selection of likely cluster members, however, was 
based on the uncalibrated color-color diagram. 

Our original region, selected to include the giant branch, was M < 18, C — M > 0.2 and 
M — T2 > 0.5. The stars observed in 1996 with photometry from our WIYN mosaic are shown in 
Fig. [J in the calibrated C — M vs. M — T2 color-color diagram. We have adopted E { B — V ) = 0.02 
(|Harrislll996l ) for M92 and have taken the reddening ratios from lHarris Canternal (|1979l ). These 
stars were selected on the basis of their position in the color-magnitude diagram. The open symbols 
are the non-members and the filled symbols are the spectroscopically-determined members. We have 
also inclu ded the loch of the [Fe/H] = through [Fe/H] = —3 stars, in 1 dex steps in [Fe/H], from 
Fig. 7 of iGeisler et al.l (|199ll ). The separation by metallicity is quite clear. We based our final 
selection of stars for the 1997 observing season on their location in this diagram. The selected 
region is indicated by the dashed box. 

Figure [2] shows the stars observed in the 1997 observing season. The lines and points are as 
in Fig. [TJ These stars were primarily selected to lie outside of 3'. We observed 162 new stars 
based on their Washington colors and 152 yielded useful velocities. Of these stars, 111 (73%) 
proved to be members. Outside 3', we measured velocities for 140 stars. Of these, 99 (71%) were 
members. By contrast, for stars observed in 1996 defined the same way, (r > 3' and with subsequent 
WIYN photometry), only 59 of the 142 stars (42%) proved to be members, although they were, 
admittedly, selected on the basis of an inferior CMD. Nonetheless, the 1997 sample extended to a 
larger distance from the cluster center and would have had a lower success rate had the candidates 
been selected in the CMD alone. Figure [3] shows the M vs. M — T2 CMD with members and 
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Fig. 1. — Washington system C — M vs. M — T2 color-color diagram for the stars observed in 
1996. Only stars with the higher-precision, WIYN photometry are shown. The filled symbols are 
the stars which were subsequently determined to be members of M92 and the open symbols are 
non-members. The solid lines are locii of common metallicity from Geisler (1991). Based on this 
diagram, we selected the region within the dashed lines for our 1997 observations. 
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0.8 1 1.2 1.4 



(M-T 2 ) 



Fig. 2. — As Figured] for the stars observed in 1997. This stars in this sample were selected only 
on the basis of their positions in this diagram. 



-8- 




0.5 1 1.5 2 

(M-T 2 ) 

Fig. 3. — Washington system color- magnitude diagram for all stars with measured velocities. The 
filled symbols are confirmed members, with the circles for the velocity variables and the squares 
non- variables. The open squares are non-members and the stars doubtful members. Most of the 
stars brighter than M = 15.5 were observed in 1996; those fainter are from the color-selected 1997 
sample. 
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non- members distinguished. Most of the stars in the 1997 sample (Fig. [2]) have M > 15.7 and 
M — T2 < 1.2 in this CMD, while most of the stars with M < 15.5 are from the 1996 sample 
(Fig. [1]). A comparison between the two samples makes it clear that member selection is easier in 
the Washington color-color diagram. 



2.2. Astrometry 



We calculated the astrometric positions of our stars using ve rsion 1.2 of the HST Guide Star 
Catalog (GSC) as our astrometric system (see lLasker et al.lll990l ). About 9000 stars in 45' x 45' 
region around M92 from the 'Quick V" Palomar Schmidt survey were used as secondary astrometric 
standards. Plate solutions were calculated for each of the DDO 51 WIYN frames in the mosaic; 
the narrow-band filter avoids any problems with differential refraction. The astrometric solution 
was good to better than the 0'.'3 required for Hydra. 



2.3. Velocity Observations 

We obtained all observations with the same spectrograph configuration we used for M15 



(jPaper I ). We used the echelle grating with an order centered at 515 nm, in the neighborhood 
of the Mg b lines. Approximately 20 nm of the order was imaged on the 2048 pixel long CCD 
for a dispersion of about 0.01 nm/pixel. The comparison source was a Th-Ar lamp. The Hydra 
multi-fiber positioner and WIYN bench spectrograph were used. The positions for the sky fibers 
were determined from the Quick V frame. 

We observed M92 over the course of four observing runs in 1996 May, 1996 June, 1997 June, 
and 1998 June. Table [2] contains a log of the observations including the number of stars observed 
in each config uration. Th e 1996 May runs observed fewer stars than the later runs since the proper 



motion list of iReed (|1992l ) was confined to the inner 13.3' of the cluster. The Schmidt images and 
the WIYN data covered a much larger area, allowing us to observe more stars per configuration. 
As the 1997 June run progressed, we did a rough reduction of each night's data and removed stars 
showing grossly different velocities from that of the cluster. This allowed us to concentrate on the 
stars more likely to be members. Most of the removed stars were in the outer part of the cluster, so, 
since it is harder to place fibers in the more crowded inner region, the number of stars observed per 
configuration went down as each run progressed. The proportion of members increased, however. 
In 1998 June we observed one configuration in order to refine our velocities of a number of stars, 
but more particularly to obtain additional spectra of a few stars with ambiguous membership. 

In all we obtained 1204 useful spectra of 596 stars. Of these 306 turned out to be definite 
members as discussed below. We include in this total velocities for six faint stars for which we 
could measure only a single velocity by combining all of their spectra. Five of these turned out to 
be members. For two of these members an additional, confirming, velocity was obtained in 1998 
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June. 



2.4. Velocity measurement 

We observed three stars, IV-10, 111-13 and 220 repeatedly. IV- 10 was our primary standard 
and was stable, as was 220, but 111-13 showed the classic jitter of bright giants. We used IV-10 as 
our velocity template for the velocity determination by cross-correlation. 

The data were reduced using the dohydra reduction package in irafH. Each observation was 
accompanied by one or more 5-minute exposures of an incandescent lamp (a "flat") taken with 
the fibers in the same configuration as the observations. Generally, configurations observed at the 
ends of the night had multiple flat exposures, but, because of the overhead involved with flats 
and especially with reconfiguring the fibers, we usually only took a single flat exposure. There 
did not appear to be any disadvantage to using single flat exposures since cosmic rays were not 
a great problem. The program exposures and bracketing Th-Ar lamp exposures were extracted 
and then divided by the extracted lamps. In the two 1996 runs and the 1998 observation no sky 
subtraction was required since they were observed in dark time. Our 1997 June run took place in 
brighter conditions and some sky subtraction was required. The sky contribution was small and 
only affected the observations at the beginning of each night. 

Cosmic ray removal was done through the simple expedient of using the IRAF continuum task 
to replace with the continuum fit all pixels more than 4 standard deviations above the fit. Unlike 



m 



Paper J . we did the cosmic ray removal before the dispersion correction and resampling. This 



prevented the negative spikes we saw in the M15 data. 

As in our M15 observations, the wavelength calibration was done using 34 to 36 comparison 
lines. The fifth-order dispersion solution generally had RMS residuals of less than 10 -4 nm or 0.05 
km s _1 . During dispersion correction the spectra were re-sampled into 2048 logarithmically spaced 
bins covering 20.7 nm in total. 

In 1996 May, 1997 June, and 1998 June we took multiple exposures in each configuration. 
The resulting spectra for each star from each configuration were summed together to produce the 
fi nal spectr a for cross-correlation. The velocities were calculated using the fxcor task in IRAF. As 



m 



Paper II . we excluded all spectra with cross-correlation peaks smaller than 0.2. These had low 
signal-to- noise ratios and gave highly discordant velocities in most cases. The template consisted 
of the sum of all exposures of star IV-10 in 1996 and 1997 suitably shifted to a common heliocentric 
velocity. This spectrum is the result of 28 hours of total integration time. 



Our error analysis proceeded as in iPaper 11 following the method of iPrvor. Latham. Hazen 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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(198 8|). We used the repeated observations of stars to establish the value for the constant in the 



Tonrv &: David (|1979l ) formula for the velocity uncertainty. In this instance we obtained C = 
12 .8 ± 0.5 k ms" 1 . This is quite consistent with the value of 13.1 ± 0.5 kms -1 we obtained for M15 
in iPaper J . As the experimental setup is the same, and M15 and M92 have similar metallicities, 



this is not surprising. 

The velocity zero point was established using 18 spectra of the twilight sky taken during 
various observing runs. During our 1997 June run sky observations were taken both at evening 
and morning twilight. We extracted the spectra in the usual manner and cross-correlated them 
against the template spectrum. The individual velocities were examined for fiber-to-fiber variations. 
Nothing significant was seen. Using the value of C above, the individual velocities have errors of 1 
kms -1 . This is primarily an effect of the mismatch between the low-metallicity template and the 
solar spectrum of the twilight sky. The mean velocity in each image had a standard deviation of 
only 0.1 to 0.2 kms -1 . These means were used in estimating the zero point. 

We noticed a systematic difference between the 1997 June evening and morning exposures 
of about 0.54 ± 0.02 kms -1 , excluding one anomalous value. We attribute this to the rotational 
velocity of the Earth, which, at the latitude of Kitt Peak and for the times of the observations, 
amounts to 0.64 kms -1 . Subtracting this contribution leaves a difference smaller than the errors 
in the means for the two twilights. We show the results in Fig. |4l The upper panel shows the 
velocities corrected for all motions except the rotation of the Earth against the fraction of a day 
from mean noon at Kitt Peak. Therefore, we corrected all the sky velocities to the heliocentric 
frame. There were still some residual differences between evening and morning observations, as 
seen in the lower panel of Fig. UJ These were 0.07 kms -1 from their mean in 1997 June, equivalent 
to about 0.02 of a pixel at our dispersion. The highly consistent behavior leads us to suspect 
some remaining systematic effect we are unable to identify. We are similarly unable to account 
for the large deviations from the mean of a couple of observations. Our zero-point is taken as the 
unweighted mean of all the corrected sky observations and is —119.6 =L 0.1km s . 

Figure [5] shows the velocities for all the stars observed. The group near —120 kms -1 are 
principly cluster members; the Galactic disk stars cluster near zero. As in M15 we see a few stars 
with velocities significantly more negative than the cluster. These are presumably field halo stars. 

We established mem bership b y requiring velocity coherence and line strengths appropriate for 
low-metallicity giants. In lPaperll . the metallicity criterion was simply an estimate of the equivalent 
width of the Mg line at 518.3 nm. For our M92 observations, however, there were several stars with 
weak spectra for which this technique gave am biguous results. Instead, we have used a variation of 



the method of lRatnatunga Sz Freeman! ([1989 ). They showed that the height of the auto-correlation 
function for a spectrum is related to the star's metallicity. Since we simply desire an uncalibrated 
measure of the metallicity, we have simplified their procedure somewhat. We have filtered the 
spectra in the Fourier regime to remove the highest frequencies before autocorrelation and have 
normalized the autocorrelation height by the log of the mean number of counts in the spectrum. 
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Fig. 4. — (top) Relative velocity of twilight sky exposures corrected for all motions except rotation 
of the Earth plotted against fraction of a day from mean noon at Kitt Peak, (bottom) After 
correction for terrestrial rotation. 
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Fig. 5. — Radial velocity vs. distance from the center of the cluster for all observed stars. The 
cluster near —120 kms -1 are mostly cluster stars. A few presumable halo stars are visible at more 
negative velocities. 
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Fig. 6. — Normalized autocorrelation height, h±, plotted against velocity. The box contains the 
locus of the members. The open circles are the four stars rejected for the reasons discussed in the 
text. The filled circles are other doubtful members. 
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We plot in Fig. [6] this h± index against the radial velocity of the stars. 



We make the following cuts in Fig [6] to establish our members. First, we reject all stars with 
velocities differing from the template by more than 35 kms" 1 . Second, we reject all stars with 
hi > 0.25. This leaves us with 310 stars within the box in the figure. Of these, four stars (open 
circles in Fig ED have been rejected from consideration for the dispersion analysis. All show highly 
discrepant velocities compared with other stars at the sam e distance and all but the last have zero 
probability of membership based on lTucholke et al.l (|1996l ) proper motions. Other stars with zero 
probability have been included in the sample if they otherwise satisfy the criteria above. 

The remaining 306 stars are listed in Table El Our J2000.0 coordinates are given in columns 
(1) and (2), followed by a stellar ide ntification. Where pos sible, we have used identifications which 
are generally derived from Table 3 of lTucholke et al.1 (|1996l ) , although we have made some additions 
and corrections where w e have deemed them appropriate. These ide ntifiers primarily come from 
Sand age fc Walkerl (119 66) (hyphenated Roman numeral prefi xes or 'x') , Cudworthl (119761) C'C ' ), and 



Buonanno et al 



(119831) f'Bu 



('N'), IZinn. Newell fc Gibson 



with addition al stars from iBarnardl (|193ll ) (' B'), 



(|l972 ) ('ZNG'), ISawver Hoggl (|l973l ) ('V'), and |Reej (|1992l ) ('R 



Nassau 



(|1938l ) 



Where available, other names from these sources are given in column (11). Where no previous 
identification exists we have give our own number. The next column contains the distance from 
the cluster center, which we take as a(J2000.0) = 17 h 17 m 07!02, <5(J2000.0) = 43°8 / ll'. / 4. The next 
three columns give our Washington photometry. Column (8) gives the number of velocities for the 
star. Where there was only a single observation, the identification of the configuration in Table [2] is 
given. The two stars with configuration 'Z' have single velocities based on sums of all the available 
spectra. Columns (9) and (10) give the weighted mean velocity for each star and their uncertainties. 
Finally, notes for the individual stars are given in column (12). Table U] contains the information 
for the doubtful stars and Table [5] (available on-line only) contains a similar list of non-members. 

Twenty-six of the members are flagged as variables for various reasons. For 15, multiple 
radial velocity measurements show them to be radial- velocity variables. This is defined as having 
probabilities of less than 1% that their x 2 are consisten t with no y aria tion. Eight stars, incl uding 
two of the velocity variable s , are i den tified in the li s ts of iReed (|1992l ) or lTucholke et al.l (| 19961 ) with 
variable stars from iNassaul (|1938l ) or ISawver Hogg! (|l973h . Most of these are RR Lyr variables. 



A further group have been flagged due to the velocity "jitter" seen in the brightest giants 
(|Cunn fc Griffinlll979h . Of the 11 giants in Fig. din the clump with M < 12.8 and M-T2 > 1.55, 
8 have multiple observations and 6 of these are v elocity variables. If we include an addition 0.8 
kms -1 velocity uncertainty ( Gunn Sz Griffin! Il979l ) to allow for this jitter, then the velocities are 
consistent. We have removed these 11 stars from consideration for our velocity dispersion analysis. 

In addition to the four discrepant stars, we also include in Table[5]three further stars. The stars 
1016 and R644 have accordant velocities, but slightly high hi values, and have been thus rejected 



3 Stars identified by 'Bu III-nn' are those listed in Table III of that paper. All others are from Table II. 
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Fig. 7. — Velocties vs. radialdistance for the members of M92. Circles represent mean velocities of 
the variables. Note the existence of stars with large velocities at large radii. The line is the mean 
velocity. The lighter 'x's are the doubtful members. (The seventh object is in the clump near the 
mean velocity at 1'.) 
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from the me mbership list. On the other hand, star 1016 has an 81% probability of membe rship 
according to iTucholke et al.1 (|1996l ) and R644 a 99% probability of membership according to iRees 
()1992l ). It is probable that they really are members, but we have excluded them from the analysis 
for consistency. Finally, the star VI-7, which is the star at 7 ' with velocity —140 kms -1 has been 
rejected based on its peculiar DDO colors in lNorris Zinnl ([19771 ). 



The final result is that we have 306 stars we consider to be members and 7 stars which are 
doubtful members. We show the velocities of theses stars against their radial position in Fig. [7[ The 
open circles indicate the 26 probable velocity-variables mentioned above, while the fainter points 
are the doubtful members. The median error for these velocities is 0.35 kms -1 with 92% having 
errors less than 1 kms -1 . 



Our sample of 306 stars includes 34 of the 35 stars for which ISoderberg et al.l (|1999l ) have 
reported velocities. Comparison of our velocities with theirs indicates a very similar mean velocity 
for the stars in their sample. However, the individual velocity differences from the mean show a 
systematic deviation between the two studies, viz. their values for this difference are about 75% 
of what we find. Their template was constructed from an average of all their bright giants and 
includes strong interstellar Na D lines as well (Pilachowshi, private communication). To the extent 
that a stellar spectrum is dominated by these interstellar lines, its cross-correlation velocity will be 
pulled toward the template velocity, i.e. to the cluster mean. Thus, the measured cross-correlation 
velocity will be a weighted average of the cluster mean and the true velocity of the star. Since 
our spectral region does not contain strong interstellar lines, our velocities are likely to be more 
reliable. 



3. Bayesian Analysis 



3.1. General Principles 



As in 
principles ( 



Paper 11. we derive the velocity dispers ion profile from our velocities using Bayesian 



Sivial Il996l: I Javned 120031 ; iMacKavl 120031 ) . Our method here differs from, and improves 



upon, that used in lPaper J . The improvements are procedural and conceptual, as will be discussed 
below. 

Our approach is to consider a set of specified, parametric, models for the velocity dispersion 
profile. These do not exhaust the full range of possible profiles, but do provide a plausible starting 
set. We then infer the posterior probability distributions for the model parameters for each model 
on the basis of our observed velocities. We also calculate the relative likelihoods of the various 
classes of models with respect to one another. This allows us to identify the model best supported 
by the data. 



Generically, if the vector p represents the parameters for a model or hypothesis H, and D 
represents the available data, then, given background information /, Bayes' Theorem states that 
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the posterior probability for p is given by 



The quantity p(D\pHI) = L(p) is the likelihood which, for the present problem, we take to be 

N 



N 

L( P ) = (2^/ 2 n^ 1/2 exp 



9 X) K - v f k i 



2 

i=l 



(2) 



where, if is the uncertainty in velocity Vi and <r(rj) is the dispersion appropriate to a star at 
radius n, 

k^^ + ainf]- 1 - (3) 
The connection with maximum likelihood methods is discussed below. 

The components of the parameter vector are the mean velocity, v, and the various parameters 
required to define cr(r). The mean velocity is a nuisance parameter in terms of inferring the velocity 
dispersion profile, and can be removed by marginalizing over it. 

The prior probability for p is given by p(jp\HI). For the most part the parameters for the 
models we use are logically independent, so p(p\HI) will factor. The denominator, p{D\HI) = 
J p(D\pH I)p(p\H I)dp, gives the overall probability of the data given a hypothesis which is vari- 
ously termed the 'evidence' or the 'model likelihood'. For a particular model this is just a normal- 
ization factor, but application of Bayes' Theorem once again gives p(H\DI) oc p(D\HI)p(H\I) for 
the posterior probability of model H. When we wish to compare two models, Hi and H 2 we can do 
this by taking the ratio of their posterior probabilities: 

p(Hi\DI) = pjDlHjl) pjHtf) 

p(H 2 \DI) p(D\H 2 I) p(H 2 \lY 1 1 

The first factor on the right hand side is the likelihood ratio using the values of p(D\HI) we 
just compute d. The second factor is the ratio of the prior probabilities for the two models. Our 



procedure in iPaper II had the shortcoming in that instead of calculating p(D\HI), we compared 
the maxima in p{p\DHI). This compares only the best set of parameters for each model. The 
proper thing to consider is the relative merits of the models for all combinations of parameters. 
Comparing the probabilities of the best sets of parameters favors a model having a narrow range 
of parameters that fit the data particularly well over one that has a broader range of parameters 
that fit the data somewhat less well, whereas the opposite ought to be the case. The issue is the 
choice of the best model class, not simply the best set of parameters. 

In this we differ from maximum likelihood methods. The latter can be derived from the 
Bayesian method in the limit where the prior probabilities are constant. The objective is then 
to find the maximum, i.e. the mode, of the likelihood function and to report the function based 
on this set of parameters as "the velocity dispersion profile", in much the same way as we did in 



Paper J . The width of the likelihood function can be used to give uncertainties for the parameters, 
but this method has no way of assessing the relative merits of various possible descriptions of the 
data. 



At the same time, it is often difficult to illustrate fully the complexity of the posterior prob- 
ability function. For situations with a large number of parameters, Markov-chain Monte-Carlo 
methods are practical, but can take a long time to generate sufficient samples to get ensure that 
the posterior is being properly sampled. For problems with relatively few parameters, direct cal- 
culation of the probability is a practical proposition. Consequently, in this paper, we compute 
the joint probabilities p(p\DHI) directly on a multi-dimensional grid. The highest dimensional 
grid we used here has five dimensions, but the number depends on the model as discussed in the 
next section. For models with radial binning, the parameter space naturally factors into smaller, 
tractable, sub-spaces. In each dimension we used a grid of 50 points concentrated on the region 
of significant probability and is well resolved. We then extract the information in this function in 
several ways to highlight not only its maximum, but also certain aspects of its shape. 

To this end we calculate the marginal probabilities, p(pi\DHI), of the each of the parameters, 
Pi by integrating p(p\DHI) over all the other parameters in the problem. Numerically, we present 
our results in terms of the moments and mode of these marginal distributions. For each parameter 
we give the mean defined as the expectation value 



We also estimate the mode of the smooth marginal distributions by constructing a parabolic fit to 
the mode and the two bracketing points on the function grid. Deviations between the mode and 
the mean highlight asymmetries in p(pi\DHI). In such cases we will use the following notation 
to indicate both these values: [x] = (x) (x) ±e x , where the value in parentheses is the mode. 
Otherwise, we will just give the mean and la standard deviation as usual, although this should not 
be taken to necessarily indicate that the marginal distribution is Gaussian. 

It should be noted that the modes of the marginal distributions correspond to the overall 
mode of p(p\DHI) only if the joint probability is symmetric. In practice, there are often significant 
correlations between parameters as highlighted by covariances calculated from joint probabilities 
of pairs of parameters. In order to highlight the range of parameters permitted by p(p\DHI) we 
give two graphical representations of our results in addition to the values mentioned in the previous 
paragraph. First, we will show the marginal probability distributions for the best-fitting models. 
We also calculate p(a\rDHI) by integrating p(p\DHI) over regions of constant cr(r). From these 
we will show the mean value and uncertainty as defined by equations ([5]) and ([6]), as well as the 
mode. Again, the difference between the mode and mean highlight asymmetries in p{a\rDHI). In 
practice, these curves will also deviate from those based on the mean or modal parameters for the 




(5) 



and the standard deviation 




(6) 



-20- 



model, again because of asymmetries in the overall probability distribution. For ease of graphical 
comparison with later models, we shall tabulate some of the points on these curves. We wish 
to emphasize that the errors on these points are by no means independent. Comparisons between 
future models and these data are best undertaken with the original velocities, not with these curves. 



In light of these new, improved, procedures, we will also reanalyze the M15 data from lPaper I . 



3.2. The Models 



In accord with the discussion in iDrukier et al.l (|2003l ) highlighting the underlying difficulties 



in inferring global properties from a set of individual velocities, the velocity-dispersion profiles 
to be derived below depend on the following assumptions. We assume that the velocities can 
be perfectly described by Gaussian distributions with a common mean and a dispersion that can 
vary with radius. Any higher-order moments will be ignored in our analysis. We assume that the 
uncertainties derived above for our radial velocities are accurate and precise and that the noise 
in the measured velocities are distributed as Gaussians with zero mean and standard deviation 
equal to these uncertainties. We assume that there are no non-members present in our sample. 
We do not exclude the possibility of there being unbound stars. Indeed, the use of a Gaussian 
implicitly assumes that the extreme tail is populated by unbound stars. We also assume that all 
the stars in our final sample are single stars with the same mass. It is possible that there are 
unidentified velocity variables, including binaries, in our sample. These would serve to inflate the 
dispersion. These provisos are included in the background information I. One model which could 
be considered is one where the distribution of the velocities is consistent with the observed errors. 
Such a model, however, is strongly rejected with respect to all others, indicating that we detect the 
intrinsic variance in the cluster velocities. We now list the models we will consider in this paper. 
This list is not exclusive, nor exhaustive, but does cover likely variations^ 

Model B: The stars are sorted by distance from the cluster center and divided into bins. The 
velocity dispersion is taken to be constant in each bin. This is the traditional approach, which we 
augment with our Bayesian estimation of the velocity dispersion. The parameters are the mean 
velocity of the sample, v, the number of bins M, and the set of dispersions, {crj}, j = 1, . . . , M. The 
choice of binning leads to various sub-types of this model class. Bins can be arbitrary, or divided 
roughly equally in terms of number of stars per bin or in terms of the radial range covered. An 
alternative is to let the numbers of stars per bin, or equivalently the radial ranges, be an additional 
set of M — 1 parameters, rj, i = 1, . . . , M — 1. 



Which is not to preclude other, potentially interesting models. The iKing (1966) model, has often been used 



to 

describe globular clusters. With its energy cut off, however, it explicitly cannot account for escaping stars, which 
would then have to be accounted for separately. Such two-component models, as discussed at the end of this section, 
are beyond the scope of this paper. Sim ilarly, we do not investigate the models including non-Newtonian effects 



( Scarpa. Marconi, fc Gilmozzil (20031 . 120061 ). 
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Model P: The velocity-dispersion profile has the form of a power-law a{r) = a\r a . The 
parameters of this model are: v; a\, the velocity dispersion at r = 1'; and a. 

Model C: The velocity-dispersion profile has the form of a power-law with a core: 

{o n a/2 
1+ r/rp M 
: , -2 r • ( ? ) 

The parameters are as in model P with the addition of the scale radius ro of the core in arc minutes. 
Note that, as in model P, o\ is normalized to be the value of the dispersion at r = 1'. This removes 
the correlation between ro and a\ that otherwise would be present. Model C is identical with model 
P in the limit rg — ► 0, and so is a more general case. For large r, a oc r a . 

Model PB: This model is a hybrid with a power-law inner section and an ou ter sectio n where 



the velocity dispersion in constant. This model is motivated by our discovery in lPaper II that the 
velocity dispersion appears to increase in the outer part of M15. It seems sensible then to consider 
that these stars may not share the velocity dispersion properties of the rest of the stars in the 
cluster. The parameters here are the same as in model P, with the addition of the radius r\ 
dividing the two populations and the dispersion, <7 out , for these outer stars. 

Model CB: This model is as model PB, but for the cored power-law. It has the parameters 
of model C, plus n and cr ou t. 

In both models PB and CB the assumption is that there are distinct radial regions containing 
distinct populations of stars. The inner population are the stars bound to the cluster, while the 
outer population consists of those stars which deviate from the cluster distribution because of tidal 
perturbations, ejection from the cluster core, or other, unknown reasons. A more realistic model 
would be one where the two populations are intermingled with additional parameters, including the 
membership probabilities for the two populations, describing the composite distribution. Similarly, 
the possibility of undiscovered binaries could be included. These are a more complicated proposi- 
tions to implement than the models presented here, but may be r equired to ful ly understand the 



data. A description of a similar approach can be found in §21.5 of iJavned (|2003l ). 



3.3. Priors 



We now consider the prior probabilities for the various parameters. The mean velocity is 
common to all the models and we take p(v\HI) to be uniform. Its normalization will drop out of 
all model comparisons, so we can safely use a flat, improp er prior. For the various dispersions, the 
correct prior is the Jeffreys prior (|Sivialll996l ; |Javnesll2003l ) p(a\HI) oc a -1 . This is equivalent to a 
flat prior in log(er). We will need to set sensible bounds, <7_ < a < <7+, on this in order to make it 
normalizable. We have used cr_ = 0.3 kms -1 , since this is the median error on our velocities, and 
smaller dispersions are likely to be meaningless. For the upper limit we use values appropriate to 
each cluster. For the models using the power-law index a, we have taken p{a\HI) to be uniform 
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over —2.5 < a < since the velocity dispersion should be a decreasing function of radius, but not 
too strongly. The lower limit is set by consideration of the range of physically sensible values for a. 

For models C and CB, since we interpret ro as a scale parameter, the correct prior is the 
Jeffreys prior. It is meaningless to speak of a core radius outside the radial range of the data, so 
the limits on ro are taken with this in mind. 

Finally, for models PB and CB, as well as variants of model B with varying bin sizes, we 
take p(ri\HI) to be uniform over the appropriate radial range. In this study we will consider two 
different sorts of binnings. The first is a standard arrangement, with the stars divided into bins of 
equal number. The second is to set a minimum number of stars per bin and then scan through all 
possible combinations of break radii in order to find the optimal bin size. We will denote models 
with these two binnings as £?jv=m and B$=m respectively, where M is the number of bins. A 
scanning approach similar to the B$ models is taken to the break radius r\ in models PB and CB, 
and, again, the prior is uniform over a specified radial range. 

Note that for models in class B with M > 2, the prior probabilities on the positions of the 
bin boundaries are not independent, since they are sequential and we require each bin to contain 
a minimum number of stars. For M = 3, for example, p(r\r2\B s=sl) = p{r\\Bs = ^I)p{r2\riBs=^I) 
follows from the product rule of probabilities, with similar expressions for additional bin boundaries. 

The limits on the various parameters for the M92 and M15 data sets are given in Table [6J 

The selection of priors for the models themselves is a somewhat more difficult proposition. In 
some cases, we may indeed have pertinent prior information. For M15, for example, the surface- 
density and surface- brightness profiles are power-law in form to within 2" of t he center of the cluste r 



(jSosin &; King]|l997l ). Inspection of the velocity-dispersion profile in Fig. 9 of lGerssen et al.l (|2002l ). 
suggests that the break in the slope of the velocity dispersion occurs somewhat further out, at 
around 1' to within a factor of 2. This difference is to be expecte d. If, for example, we con sider 



the single-mass Wq = 9 King (1966) model shown in Fig. 4-11 of iBinnev Sz Tremaind (| 198 71 ). we 
can see that the line-of-sight velocity dispersion drops to half its central value at about 30r c . We 
could consider this an effective core radius for the velocity dispersion profile, r c ,v Thus, it would 
not be unreasonable to expect to see some evidence for flattening within our data sample. If the 
core radius in the velocity dispersion is unresolved, it may be difficult to differentiate between a 
true core, and the situation we will see in Example PL+F of §3.51 below. Taking this into account, 
p(P\Imi5) /p(C\Imw) ~ 5 would not be unreasonable for the two options of a power-law (model P) 
or a power-law with core (model C). This expresses our relative ambivalence with respect to the 
two models, but a desire for some stronger evidence of flattening before we are prepared to believe 
it. Like any prior, it is a statement of how strongly we hold any given proposition. 

For M92, our preferenc e is for model C over model P, sin ce the cluster shows a clear core radius 



of 14" in the density data (jTrager. King. &: Djorgovskilll995l ). We would then expect r c ,v ~ 7'. In 



terms of the parameterization of model C, ro = r c>v /yA l l a — 1. For a < —0.25, ro > 1', and 
is well within our data sample. Thus we expect to see clear evidence for a core in the velocity 
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data. In principle, this argument can lead to a stronger joint prior on ro and a than the relatively 
uninformative priors proposed above, but we shall refrain from doing so here. Instead, we take 
p{C\lM92) /p{P\Im92) = 100, expressing our strong prior expectation of model C. 

In the absence of suitable prior information, the same considerations which lead to the nu- 



merical values above suggest we employ the "device of imaginary results" of iGoodl (|l950l ). Before 
doing the analysis of the results, we ask ourselves how much better would the likelihood of one 
model have to be over another for us to consider the two models equally probable. The inverse of 
this ratio is then the desired prior. If, for example, the possible models for a set of globular cluster 
data are either one with or one without a central black hole of unknown mass, and, being skeptical 
of the proposition, we would need the likelihood in support of a black hole to be 50 times that 
without, then the proper prior ratio is p{B\DI) /p(B\DI) = 1/50 where B is the proposition that 
the cluster contains a black hole and B its negation. 

The case of binned models is interesting in this regard. No one actually thinks of a globular 
cluster as being stratified in the way implied by binning. The main reason to use such a model is 
to look for features in the profile not described by the simpler parametric models such as P, C, 
etc. Increasing the number of bins allows for finer detail in the resulting profile, but also decreases 
the precision of the measurement of the dispersions. The Bayesian formalism naturally takes this 
effect into account. A model with more bins is penalized by the additional priors attached to 
additional parameters. If, nonetheless, the total probability is maximized for a particular set of 
bins, this is a way to optimize the choice. Such optimization is, of course, subject to the prior on 
the number of bins, but we shall take that to be uniform up to some maximum number of bins. Our 
objective with this class of models is to look for the most probable binning in case there is some 
obvious feature not captured by the other models under consideration. Nonetheless, given the lack 
of physical motivation for a binned model, we set p{B\I) to be 0.02 times the prior probability of 
the preferred model of P or C as appropriate for M92 and M15. This probability is intended to 
apply to the entire collection of binned models reflecting the flexibility the observer has in choosing 
the binning. In the subset of models with variable-width bins we require a minimum of 20 stars 
per bin for M = 2, 3 or 4. (The calculation becomes very long for higher values of M in this case.) 
Applying the same limit to bins with an equal number of stars per bin, permits all models with 
M < 14 for the M92 data and M < 10 for the M15 data set. Dividing p(B\I) equally between the 
17 cases for M92 and 13 cases for M15 gives relative prior probabilities of 1.2 x 10 -3 and 1.5 x 10~ 3 
per binned model for the two clusters. 

Finally, we consider the odds to be even on whether the outer zones in models PB and CB 
exist or not. 

At this point, we invite the reader to decide what prior ratios are appropriate for the various 
models we describe in £|3.2l in the cases of M92 and M15. Below, we will present the results in such 
a way that the posterior odds can be recalculated for other choices of prior ratios. 
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3.4. Implementation Details 

The posterior probability distribution for each model was calculated on an ^-dimensional grid, 
where N is the number of parameters in the model. The marginal probability for each parameter 
was calculated from the values on this grid by numerical integration over all the other parameters. 
Along the way, the two-dimensional marginalizations were examined and used to calculate the 
covariances between pairs of parameters. To increase the efficiency of the calculations, we first 
conducted a pilot study using a coarse grid covering the full ranges of the priors. This was used to 
locate the region of non-trivial probability in each parameter, which was then covered by a finer grid 
in our final calculation. The marginal probability distributions were saved. We also calculated the 
first and second moments of these to give the expectation value and variance of the distributions. 
In cases where the posterior probability distribution is close to a Gaussian, these are close to the 
mean and variance of that Gaussian. 

In the cases where one or more break radii are themselves parameters, the above calculation 
was done for each set of break radii and the total probability is naturally the marginal probability of 
that set of radii. The posterior probability distributions of the break radii are highly non-Gaussian 
and rather than specifying the rather unenlightening moments, we use the mode as the typical 
value. The other parameters were marginalized over the break radii. 

3.5. Examples 

We now proceed to use some mock data sets to demonstrate the use of our methods. In all 
cases, the radial distributions and velocity uncertainties are those of our M92 data. For each mock 
data sets we used our Bayesian procedures to estimate the parameters for all the models discussed. 
For the binned models we used equal number bins running from 1 to 20 bins and scanned bins 
with 2 to 6 bins, i.e. with 1 to 5 breakpoints, and a minimum bin size of 40 stars. We present our 
results in Table [7J The log of the likelihood ratio with respect to the best model is listed for each 
of the mock data sets. For each of the binning schemes, the total probability presented is the sum 
of the total probabilities for all binnings. We also give the results for the best value of M in each 
binning scheme. In all cases, the mean velocity of the mock data sets was consistent with the input 
value. Although we have dropped the units for simplicity, they should be clear from the context. 
We discuss each mock data set in turn. 

PL: A power-law with o\ = 6.5 and a = —1. The best model is P with (<7i) = 6.0 ± 0.4 and 
(a) = —1.05 ± 0.06. Model C is 50 times worse, and, since the best core radius is at the lower limit 
of the prior, this supports the primacy of model P in this case. The only comparable model is PB, 
where the o\ and a probability distributions are nearly identical. In the outer region, p(a ont \DI) 
is strongly peaked at the lower limit of the allowed range. In other words, the results for model 
PB are effectively identical with those of model P, but the model is penalized by its additional 
parameters. The steep gradient makes the constant dispersion within each of the bins in model B 
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a very poor representation of the data. 

CPL: A cored power-law with with o\ = 6.5 and a = — 1 and r$ = 1. The best model is C 
with (a x ) = 6.2 ± 0.5, (a) = -1.0 ± 0.1, and (r ) = 1.1 ± 0.4. The next best model is CB, but the 
situation here is the same as in the previous case. The parameters in the inner region are the same 
and the fit is penalized by the existence of the additional parameters. 

PL+F: A power-law with o\ = 6.5 and a = — 1 inside a radius of 8, but beyond this the 
dispersion is flat at 2. The most probable model is CB, with PB running a very close second. 
The core radius probability distribution is strongly peaked at the lower limit of the permissible 
range (0.3), giving only an upper limit. This indicates some flattening of the dispersion profile 
towards small radii, but the core is not resolved. This is easily understood in terms of the spa- 
tial distribution of the stars. In the innermost sampled region, only a small fraction of the stars 
have measured velocities. The probability of any of these stars having a high velocity is small, 
so the velocity dispersion appears to increase less quickly at the inner edge of the data. Thus 
we can see hints of a core radius, even when the model underlying the mock data contains no 
such flattening. The estimated values of oi and a are similar in the two models, but are less well 
defined in the more general model CB as seen in Figure [HJ As is also shown, the two models 
do agree on the dispersion in the outer region. The asymmetry in the probability density distri- 
bution for <T out is typical of dispersions measured from a small number of stars. For model CB 
(<ri) = 6.4 ± 0.5, [a] = -1.2 (-1.1) ± 0.1, [r Q ] = 0.6 (0.3) ± 0.3, and [cr out ] = 1.9 (1.6) ± 0.5, while 
for model PB = 6.0 ± 0.4, (a) = -0.97 ± 0.06, and [<7 out ] = 1.9 (1.7) ± 0.5. There is a minor 
disagreement as to the location of the boundary between the two regions depending on which side 
to place the three points lying between 8'.3 and 10'. 

D: A flat distribution with a = 6.5 except between 1'.5 and 3f5, where the dispersion is double. 
This is admittedly unrealistic, but is designed to show the flexibility of the procedure. The best 
model is a binned one with two adjustable boundaries giving three bins. The modal boundaries lie at 
lf53 and 3'56 and the velocity dispersions in the three bins are (cti) = 7.2 ± 0.7, (02) = 13.2 ± 1.0, 
and ((73) = 7.4 ± 0.5, with the bins containing 65, 91, and 118 stars in order of radius. The next 
best model is one with three adjustable bin boundaries, but the final result is degenerate in the 
sense that the outer two boundaries have the same mode. This model is 7 times less probable than 
the best model. The other models are at best 10 times less probable, and must be rejected for any 
sensible prior. 

S: The velocity dispersion is given by a(r) = 5. + 2. sin(r-7r/14). This mock data set is designed 
to show what happens when none of the models is suitable to the data. In this case, the most 
probable model is simply a flat distribution with (a\) = 6.1 ± 0.3. Nearly as likely is a two-zone 
model with a modal bin boundary at 6'. In this model (a\) = 5.6 ± 0.4 and [02] = 7.3 (6.5) ± 0.9. 
As shown in Figure [9] the probability distributions for this latter model are quite non-Gaussian, 
with the probability distributions for both dispersions showing hints of bimodality. Results like 
these for real data would indicate a great cause for concern. 



-26- 




Fig. 8. — Probability density distributions for the velocity dispersions (top) and power-law index 
(bottom) for the mock data set PL+F. Solid lines: model CB. Dashed lines: model PB. The 
thick vertical lines are at the input values of the three parameters. 
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Fig. 9. — Probability density distributions for the velocity dispersion for mock data set S. Solid 
line: Single dispersion model. Dashed lines: Two-dispersion model. The line at top represents the 
range of input dispersions for this data set. 
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We now turn to our velocities in M92 and M15 and analyze them in the same way. 

4. Kinematic Results 
4.1. M92 

We give the results of our Bayesian analysis of the new M92 data in Table The likelihood 
of our five model classes relative to model CB is given in the second column of the table. With the 
priors in the third column — based on the discussion in §3.31 — the posterior probability of the model 
classes is given in the fourth column. Model CB is the most likely, but model C is only a factor of 
2 less probable, so we shall consider them together. Models P and PB are an order of magnitude 
less likely than mode CB, confirming the existence of a core, and with our prior information, they 
are up to 10 3 times less probable. 

Of the -B-class models, the model with two equal sized bins is significantly more likely than 
any of the others. It is interesting to note that the model with two adjustable bins comes up with 
much the same result, shuffling eight stars from the inner to the outer bin, at the cost of a lower 
likelihood due to the prior on the position of the bin boundary. The B^ = 2 model has dispersions 
of (cr x ) = 5.8 ± 0.4 kms -1 and (cr 2 ) = 3.5 ± 0.2 kms -1 . At the bottom of Table M we give the 
likelihoods relative to model CB for the most likely binned models. Increasing the number of bins 
does not improve the match to the data. Given our prior, we reject the binned models, individually 
and as a class. 

The probability density distributions for the various parameters in models CB and C are 
shown in Figure [TU1 The probability distributions for the mean velocity and u\ are consistent 
for both models and give (v) = —121.2 ± 0.3 kms -1 and (a±) = 6.3 ± 0.5 kms -1 for model C and 
(v) = —121.1 ± 0.3 kms -1 and (a\) = 6.2 ± 0.5 kms -1 for model CB. They do differ in the pa- 
rameters giving the shape of the velocity-dispersion profile. For model CB, [a] = —0.8 (—0.5) ± 0.4 
and [ro] = 2.5 (1.5) =L 1.6 arc minutes. As can be seen in Figure [10} the probability distributions for 
these quantities are broad and asymmetric. In addition, they are anti-correlated; shallow profiles 
with narrow cores work equally well as steep profiles with broad cores. The modal value for n is 
at 10', and it appears that the stars inside this radius can only weakly constrain the shape of the 
curve. This can be seen graphically in Figure [Til (Fiducial points for the displayed M92 profiles 
are given in TableEl) The upper shaded region shows the 1-a spread in p{a\rDCBIu92)- Inside 8', 
a range of velocity profiles are possible. Outside this radius, the profile is affected by the possible 
break radius and independent velocity dispersion. (This can be seen in the step-like increase in the 
outer part of the model CB. For such models, one of the parameters is the boundary r\ and this 
parameter has a fairly wide range of values with non-trivial probability. Consequently this part of 
the profile can be either inside or outside t\ and the curve reflects this ambiguity.) For these data 
and model CB, [<r ut] = 4.0 (3.4) ± 1.3 kms -1 . The asymmetry in the probability distribution (see 
Figure [T0b). is typical of cases where there are few stars available to constrain the value of the 
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Fig. 10. — Probability density distributions for models CB (solid) and C (dotted) of the full M92 
data set. (a) Mean velocity, (b) Velocity dispersion at 1'. The dashed line shows the distribution 
for <T out in the CB model. The long tail to higher velocity is typical of dispersion distribtutions 
based on only a few stars, (c) Power-law slope, (d) Core radius. 
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Fig. 11. — Velocity disperion profiles for our M92 data plotted with the velocities. The solid curves 
are the mean of p(o\rDI) while the shading indicates the la errors based on this distribution. The 
dashed curves indicate the modal results as discussed in the text. The upper curves are for model 
CB, while the lower curves are for model C . The mean curves for each distribution are reflected 
about zero and plotted as dotted curves to ease comparison. The two stars discussed in the text 
are labeled. 
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dispersion. 

For model C, (a) = —0.4 ±0.1 and [ro] = 1.3 (0.3) ±0.8 arc min. The former is well con- 
strained, but the latter has its mode at the lower limit of its permitted range. This indicates 
that there is some flattening to the profile, but the core appears smaller than would be expected. 
The profile is shown as the lower shaded region in Figure HU Model C is strongly constrained by 
the three outer stars with large deviations from the mean velocity. In order to account for these 
velocities, a relatively shallow, yet coreless profile is favored. Model CB puts these stars in the 
outer zone, removing their influence on the parameters of the inner profile, but the velocities of the 
remaining stars can be fit by a large range of profiles. Thus, model CB is preferred. 

The most deviant velocity is that of star 1891 in our list, lying at 10^7. Under some circum- 
stances one might be tempted to reject it outright from the cluster sample. This happens to be 
the westernmost star in our sample. It is fairly faint, but lies on the red giant branch in the CMD. 
It has three r elatively low quality y elocity measurements, but these agree within their respective 
uncertainties. iTucholke et al.l (119961) give the sta r only an 8% probability of membership, but this 



is not the only star with a low ITucholke et al.l (|199d ) probability that satisfies our criteria for 
membership. A multi-component model such as that described at the end of £13.21 would directly 
incorporate borderline cases such as this by assigning some probability to their membership. For 
the present, we will just go to the other extreme and redo the analysis without star 1891. 

The results for v and o\ are the same with or without this star. For model CB, [<7 ou t] = 3.9 (2.8) ± 1.7 
kms -1 , the large skew and uncertainty reflecting the fact that this distribution is derived from only 
the two outermost stars. For both models CB and C, the probability distributions for a and ro 
are somewhat broader and shifted with respect to the distributions in Figure [TUJ For model CB 
[a] = —0.9 (—0.6) ± 0.4 and [ro] = 2.8 (1.8) ± 1.6 arc min. As in the previous case, o\ and a are 
strongly correlated with ro. For model C, (o"i) = 6.3 ± 0.5 kms -1 and (a) = —0.5 ± 0.2. There is 
now a peak to the ro distribution and [ro] = 1.5 (1.1) ± 0.9 arc min. 

From a comparison of the results of model C probability distributions for all possible subsam- 
ples obtainable by dropping a single star, that obtained by dropping star 1891 stands out. Dropping 
this star gives a significantly larger increase in the model probability than is the case for any other 
star. The only comparable subsample is the one dropping the outermost star (#2349). In only 
these two cases do the model parameters differ from those given by the full sample. And it is only 
in these two cases that the mode of the ro is not at the lower limit of the permitted range. Given 
these results, it is possible that it is only these two stars that are anomalous, rather than all the 
stars beyond a certain radius. Such a suggestion constitutes a new model class, one where a certain 
fraction of the stars are excluded, but we will not pursue this model here. We will note, however, 
that a C-class model without either of these stars gives (v) = 121.1 ± 0.2 kms -1 , (<7i) = 6.3 ± 0.5 
kms -1 , (a) = —0.7 ± 0.3, and [ro] = 2.2 (1.7) ± 1.2 arc min. From a consideration of the various 
samples and models discussed above, it appears that o\ = 6.3 ± 0.5 kms -1 , a = —0.6 ± 0.3, and 
ro = 2' ± V is a reasonable description of the M92 velocity dispersion profile. For these values of 
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a and ro we can calculate r Cjl) = 6', a value consistent with our expectations based on the density 
profile as discussed in £ 13.31 

There is little evidence from these data for a dominant population of stars in the outskirts of 
M92 with higher-than-expected velocities. We have, however, identified two stars which are highly 
inconsistent with the velocity distribution of the remaining stars even though their properties are 
otherwise consistent with cluster membership. We will discuss the possible origin of these stars 
following a reevaluation of the M15 velocities. 



4- 1.1. Rotation 



Although this paper is focused on the ve locity dispersion profile, we would like to briefly 
follow up on the claim of iLupton et al.1 (|1985l ) that the cluster is rotating. To that end, we've 
taken our preferred M92 model from the previous section, assumed it to be correct, and looked for 
sinusoidal rotational signals for six radial bins. We then subtracted off this signal and estimated 
new parameters for the profile model, iterating in this fashion until the values of the parameters 
stopped changing. The parameters for the velocity dispersion profile are consistent with those given 
above, except for a\ which drops to 6.0 ± 0.5km s -1 . The rotation profile is shown in Fig. [12] and is 
compared to that of ILupton et al.l (|1985l ). (We've used their value of r c = 0.74 pc and a distance 
of 8.2 kpc to convert their abscissae to arc minutes.) A rotational signal is only seen in half the 
bins; for the others a la upper limit is shown. The position angle is consistent for the bins where 
a signal is seen. 



4.2. M15 



In iPaper II we concluded that there was an indication that the velocity dispersion of M15 
increased at large radii. In Table [10] we present our likelihoods and posterior probabilities in the 
same way as in Table [8] In this case, the most likely model is CB, but once our prior ratios are 
taken into account, the most probable models are P and PB by a small margin. The probability 
densities are plotted in Figure [13] None of the binned profiles is at all probable. 

The higher likelihood of model CB with respect to models PB and P lies in its ability to give 
a slightly flatter profile in the center. In this case, the ro distribution is strongly peaked at its 
lower limit, indicating that whatever core there is in the velocity profile is unresolved. This is not 
unexpected, since as we discussed in §3.31 the velocity dispersion profile in M15 flattens at around 
1'. For model CB, (eri) = 8.1 ±0.8 kms -1 , (a) = -0.7 ±0.2, [r ] = 1.1 (0.3) ±0.7 arc min, and 
(cr out ) = 3.7 ±0.7 kms' 1 . Unlike M92, a and ro are uncorrelated here. The modal value for n is 



The profiles for models P and PB are shown in Figured!] The two models differ in the way they 
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Fig. 12. — Rotation amplitude for M92 as a function of radius. The first, fo urth, and fifth point s 
are la upper limits. The open points are the bins with detected rotation from lLupton et al.l ([1985). 
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Fig. 13. — Probability density distributions for models P (dotted), PB (solid), and CB (dashed) of 
the M15 data set. (a) Mean velocity, (b) Velocity dispersion at l'(o"i) and in the outer region (cr ou t)- 
The outer region contains more stars than in the comparable case in M92, so the distributions of 
(Tout are less skewed, (c) Power-law slope, (d) Core radius for model CB. 
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Fig. 14. — As Figure [TT] for the M15 data from lPaper t . The upper curve is for model PB and the 
lower for model P. 
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accommodate the stars beyond 9'. The slope and scale of the power-law are mainly set by the high- 
velocity stars of the inner cluster. For model PB (cxi) = 8.0 ± 0.8 kms" 1 , (a) = -0.47 ± 0.08, and 
(°out) = 3.8 ± 0.9 kms -1 , with the modal value for r\ at 9'. Model P is constrained to a slightly 
shallower drop ((a) = —0.39 ± 0.06) and a smaller central amplitude {{(T\) = 7.5 ± 0.6 kms -1 ). 
Nonetheless, the velocity dispersion profile at large radii for model P is less than that of model 
PB. Model PB fits the data somewhat better, but the improvement is balanced by the need for 
two extra parameters and the result is that the two descriptions remain equally probable. On the 
one hand, the data are compatible with a single power-law, but on the other, they do not exclude 
at any level the possibility that some fraction of the stars outside about 9' have a somewhat higher 
velocity dispersion than that expected by extrapolating the most probable power-law for the stars 
inside that radius. 



5. Comparison with N-body models 



In iPaper II we suggested that the apparent flattening or uptu rn in the velocity di spersion of 
M15 at large radii was due to the effect of the Galactic tidal field. iDrukier et al.l (Il999i ). however, 
showed that during core collapse, a globular cluster ejects stars from the core. It is therefore 
worthwhile to consider the effects of unbound stars on the estimation of the velocity-dispersion 
profile. While the seemingly systematic effect in M15 can be explained either way, the apparent 
high-velocity members seen in M92 might more naturally be explained by the ejection mechanism. 

We have been carrying out a separate program to simulate the dynamic evolution of star clus- 
ters using GRAPE N-body supercomputer^ at Indiana University. We have run series of isolated 
cluster models with identical stars to well beyond core collapse, as a benchmark for comparison 
with sim ulations that includ e tidal effects. Unlike our previous Fokker-Planck models of isolated 
clusters (jDrukier et al.lll999i ). full N-body models allow th e possibility of s t ars a cquiring positive 
energies, i.e. exceeding the escape velocity. As noted by Ijohnston et al.1 (|1999l ). the increasing 
dominance of the halo velocity distribution by escaping stars can produce a flattening of the halo 
velocity-dispersion profile. We examine one such model here to illustrate the effects on the velocity 
dispersion profile of unbound stars expelled from the core. 

For comparison with the M92 results, we used N-body data from a GRAPE-4 run with N = 
8192 identical point-mass sta rs. The initial state was a Plummer model with no primordial binaries. 
We used the NBODY4 code (lAarsethlll999l ) to evolve the model through core collapse and well into 
the post-collapse phase. Stellar escape occurred primarily as a result of cumulative energy increase 
from single-single scattering in the contracting core. This produces an isotropic stream of escaping 
stars. The escape rate increased as the model approached core collapse, with about 2.5% of the 
cluster mass lost to escape at the time of core collapse. As a compromise between M92 and M15, 



See lMakino et al.1 (|1997f ) for a description of the GRAPE hardware development program. 
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we selected a data snapshot from near the time of core collapse for analysis. 

In order to compare this N-body model to the M92 data, we rescale the N-body positions to 
give the model roughly the same limiting radius as M92. The velocities are scaled arbitrarily to 
give velocities with numerical values similar to the M92 ones. We then selected stars from these 
projected lists to have roughly the same radial distribution as the M92 data. We did this in two 
ways. First we selected our sample from only the bound stars. Then we reselected the sample in the 
same way, but this time allowed unbound stars to enter each sample. Each of these unbound stars 
effectively replaced one bound star in the bound sample, so the total number of stars is maintained. 
In each case a single, obvious, strongly unbound "outlier" has been rejected. As in £ 13,51 the units 
for the parameters should be obvious from the context. 

We show the results for a typical projection in Figure HSJ Stars common to both samples are 
shown as hexagrams. Filled triangles are unbound stars in the second sample. They replace the 
bound stars from the first sample shown as open triangles. The shaded regions and curves are as 
for the Figure PTT1 The best model class for the bound samples is shown at positive velocities, while 
that for the unbound samples is shown at negative velocities. The dotted lines show the reflections 
of the mean curves about zero velocity for the purposes of comparison. 

For the bound sample, models P and C are equally likely, but the ro distribution for the latter 
is peaked towards the center, supporting the slightly higher probability of model P. For model P, 
(<7i) = 4.1 ± 0.3 and {a) = —0.56 ± 0.06. We show model P. For the unbound sample, there are 
17 unbound stars out of the remaining 279. The most likely model is P>s=2 with (<ji) = 4.5 ± 0.4 
inside 1.76 and (erg) = 2.2 ± 0.1 outside that radius. A factor of three less likely, but more probable 
on physical grounds, is model PB (shown in the figure), with the outermost nine stars in the outer 
zone for the modal value of r% = 10. Of these, five are indeed unbound. The outer dispersion is 
[c"out] = 2.7 (2.3) ± 0.8, but the difference between the mean and the mode indicates a strong skew 
to higher values since there are few data points to constrain it. In the inner region, (ai) = 4.1 ± 0.3 
and (a) = —0.44 ± 0.06. The central scale is the same as for the bound sample, but the slope of the 
power-law is now somewhat flatter in order to take into account the velocities of the unbound stars at 
intermediate radii. Model P is a factor of 9 less likely with (o"i) = 3.9 ± 0.3 and {a) = —0.36 ± 0.05. 
The unbound star at a radius of 12 affects the model in much the same way as star 1891 in M92. 
Removing this star provides the greatest improvement to model P and gives (<ti) = 4.1 ± 0.3 and 
(a) = —0.43 ± 0.05, nearly identical with the previous model PB where this star is beyond r\. The 
effect of the other unbound stars remains, however. 

From this comparison, it appears that unbound stars can indeed modify the inferred velocity 
dispersion profile, generally giving a shallower slope than we find when only bound stars are used. 
Extreme outliers at large radii can be identified by the effect their removal from the sample has on 
the parameters and overall model probability. On the basis of this admittedly limited examination 
it would appear the two rejected stars in M92 are likely to be unbound stars currently exiting the 
cluster. 
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Fig. 15. — The stellar samples and most probable profiles for the iV-body model. The stars 
indicate stars selected for both the bound and unbound samples. The filled triangles are unbound 
stars which replace the open triangles in the unbound sample. The curves are as in Figure [TTJ The 
upper shading is for model P for the bound sample. The lower shading is for model PB for the 
unbound sample. 
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M15, being in deep core collapse, should have more escaping stars and should show larger 
effects, but the d etails wil l depend on the radial distribution of the unbound stars. It is possible that 
the sample from iPaper II contains some unbound stars at relatively small radii and, consequently, 
the velocity dispersion profile has been inflated in the inner region. The continuation to larger 
radii is consistent with the velocities in that region, perhaps disguising the effects of high velocity 
stars there. In any case, the ejection of unbound stars from the cluster core is a stochastic process 
and not all snapshots or projections will contain such stars, or have the same radial distribution 
of unbound stars. What this example does demonstrate, however, is that sporadic star ejection 
from the core can produce an apparent increase in the velocity dispersion in the outer region. This 
is especially so considering that, going out from the core, the unbound stars form an increasing 
fraction of the cluster-associated stars. 



6. Summary 

We have carried out an extensive investigation of the global radial velocity field in M92, 
using the WIYN telescope for both photometry and spectroscopy. We obtained high-accuracy 
(median error 0.35 kms -1 ) velocities for 299 probable cluster members, thereby greatly increasing 
the number of stars with measured radial velocities in this rich cluster, which has a very high 
density, yet resolved core. This new dat a set nicely complements that from our previous study of 



the collapsed-core cluster M15 (jPaper J ). 



We selected likely cluster candidate members for spectroscopy by several photometric methods. 
The most efficient of these involved obtaining a large-area, 3-band Washington photometry mosaic 
of M92 with WIYN. Using a photometric metallicity index, we were able to successfully identify 
cluster members with ~ 70% efficiency within the projected cluster halo (r > 3'), where the vast 
majority of stars are nonmembers. 

We've introduced an im proved B ayesian analysis scheme and have applied it to both this data 



set and the M15 data set of IPaper a . Of our models, the most probable one for the M92 data is 
one described by a core radius of about 2', velocity dispersion at 1' of 6.3 kms -1 , and an outer 
power-law slope of between —0.5 and —0.8 depending on which, if any, of two, outer, high-velocity 
stars are included. It is these stars that appear to increase the dispersion in the outer part of the 
cluster. The M15 velocities can be described equally well by both a single power-law with slope 
of —0.5 and velocity dispersion at V of 8.0 kms -1 , or by an inner power-law with slope —0.4 and 
scale of 7.5 kms -1 plus an outer region beyond 9' with a dispersion of 4 kms -1 . This region could 
be populated by some fraction of unbound stars. It is not unlikely that a few unbound stars are 
also present in the inner region of the cluster, particularly between 2' and 5'. There is also some 
evidence that the slope of the power-law flattens at the inner limit of our sample. 



Our consideration of the velocity profiles in a GRAPE-based N-body model of an isolated 
cluster suggests that our two data sets may contain a number of escaping stars that were boosted 
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to positive energies by the same internal relaxation processes that drive core collapse. Clearly, 
more data are needed to answer the questions of whether the velocity profiles do indeed flatten at 
large radii and to identify the physical mechanism behind this. For a realistic cluster in a tidal 
environment, the development of the halo must involve an interplay between internal two-body 
relaxation and external tidal influences. 

Measures of mass segregation and velocity anisotropy may provide the best means to gauge 
the relative importance of two-body relaxation and tidal effects in deter mining the struc t ure o f 
cluster halos. Mass segregation is best investigated by very deep imaging. lAndreuzzi et al.l ((2000) 
have used Hubble Space Telescope WFPC2 photometry to investigate mass segregation in M92. 
They find a significant change in the slope of the mass function between fields at 11 and 19 core 
radii (2f5 and 4f3, respectively), in the sense of a steepening with increasing radius as expected 
from two-body relaxation. It would be useful to determine the mass function, in a similar way, 
at an even greater distance from the cluster center in order to determine whether this steepening 
continues into the outer halo. Given the sharp drop of stellar surface density in the halo, this may 
require larger area, ultra-deep, ground-based imaging to obtain sufficient statistics. 

Radial velocities, alone, do not strongly constrain the anisotropy profile. Both radially-biased 
and tangentially-biased orbit distributions are capable of producing a halo flattening of the pro- 
jected velocity-dispersion profile. In c ontrast, the anis otropy profile can be directly determined 
from proper motions, either alo ne (e.g. I C ud wor t hi 1 1 9 79l . for M3) or in cqnj un ction with radial ve- 
locities (e.g. lLupton et al.lll987l . for M13). In the case of M92. ICudworthl (119791 ) finds little evidence 
for anisotropy, although it appears that this data set would only be capable of detecting a strong 
signal. 



Testa et al.l (|2000l ) have used the Digitized Palomar Observatory Sky Survey (DPOSS) to 
examine the spatial structure of M92 near and beyond the ~ 15' tidal radius. They found clear 
evidence for a flattening of the profile beyond their best fit tidal radius of rt = 12f3, which they 
interpret as a extra-tidal halo that extends to approximately 30'. The structure of this halo is 
essentially circular, showing at best weak e yidence for the sort of b i polar elonga t ion — "tails"— 



seen in some other clusters (e.g. Palomar 5: lOdenkirchen et al.l 120011 ). ITesta et al.l (12000 ) suggest 
that this may indicate that the extra-tidal halo of M92 consists of stars evaporating from the 
cluster, which have not yet formed an escaping tidal stream. This would tend to argue that two- 
body relaxation and possibly tidal shocking dominate over tidal stresses in determining the halo 
structure and mass loss rate rate in M92. 

Improved dynamical models for clusters, which take advantage of the growing datasets provided 
by high resolution imaging and spectroscopy, should provide a more sensitive measure of the role 
of tidal influences on cluster halos. This will lead to improved estimates of cluster mass loss rates 
and dissolution timescales. 
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Table 1. Comparison of M92 and M15. 



Property* 1 M92 M15 



My 


-8.20 


-9.17 


logpo (L pc~ 3 ) 


4.29 


>5.38 


rc (') 


0.23 


<0.025 b 


Th (!) 


1.09 


1.06 


n (! ) 


15.2 


21.5 


C 


1.81 


>2.5 


[Fe/H] 


-2.29 


-2.25 


R Q (kpc) 


8.2 


10.3 


Rgc (kpc) 


9.6 


10.4 


Z (kpc) 


4.7 


-4.7 



a These dat a are from th e online 
compilation of lHarrisI (|1996l ). 



ISosin fc Kind ([l997|) 



Note. - 

My = total absolute magnitude 
po = central luminosity density 
r c = core radius 

= half-mass radius 
r t = tidal radius 
c = concentration parameter 
[Fe/H] = metallicity 
Rq = solar distance 
Rgc = galacto centric distance 
Z = distance from galactic plane 
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Table 2. Observing Log for M92 



(Honfi P*n ration 


UT Date 


HJD 


"Fjxnosnrp timp 


4k 

JT 


ml 


1996 May 14 


10217.799 


2400s 


38 


m2 


1996 May 16 


10219.757 


1800s 


26 


m3 


1996 May 17 


10220.764 


1800s 


27 


m4 


1996 May 17 


10220.852 


3000s 


24 


m5 


1996 May 18 


10221.769 


3000s 


32 


m6 


1996 May 19 


10222.743 


4200s a 


29 


m7 


1996 May 19 


10222.845 


3600s 


31 


M6 


1996 June 19 


10253.685 


1800s 


29 


M7 


1996 June 19 


10253.740 


1800s 


32 


M8 


1996 June 19 


10253.785 


1800s 


31 


Jl 


1996 June 21 


10255.687 


1800s 


88 


J2 


1996 June 21 


10255.747 


1800s 


87 


J3a 


1996 June 21 


10255.818 


991s a 


85 


J3b 


1996 June 23 


10256.724 


1800s 


85 


J4 


1996 June 23 


10256.793 


1800s 


77 


J5 


1996 June 23 


10257.803 


1800s 


66 


A 


1997 June 22 


10621.724 


7200s 


71 


B 


1997 June 22 


10621.875 


7200s 


69 


C 


1997 June 23 


10622.727 


7200s 


60 


D 


1997 June 23 


10622.888 


6760s 


57 


E 


1997 June 24 


10623.708 


7200s 


57 


F 


1997 June 24 


10623.833 


7200s 


53 


G 


1997 June 25 


10624.723 


7200s 


54 


H 


1997 June 25 


10624.882 


7200s 


58 


K 


1997 June 26 


10625.725 


7200s 


55 


L 


1997 June 26 


10621.853 


7200s 


54 


M 


1998 June 16 


10980.690 


7200s 


69 



Observation affected by clouds 



Table 3. Data for cluster members 



R.A. 


Decl. 




r 










V 




Notes b and 


(J2000.0) 


(J2000.0) 


ID 





T2 


M-T2 


C-M 




(kms -1 ) 


(kins -1 ) 


Other Names 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


17 16 19.65 


+43 14 26.0 


1891 


10.65 


15.44 


1.03 


0.56 


3 


-113.08 


0.87 


T 


17 16 22.24 


+42 58 04.9 


516 


13.01 


13.62 


1.09 


0.75 


13 


-119.02 


0.09 


T 


1 i ID Z4. ID 


-r^to U / 1Z.41 


VTT Q 
V ll-o 


7 88 
1 .oo 


10.410 


J..UD 


u.oo 


9 


-1ZO. ID 


u.oz 


1 


17 16 24.38 


+43 11 45.5 


1924 


8.55 


15.47 


1.02 


0.49 


4 


-118.05 


0.44 


T 


17 16 27.68 


+43 08 55.1 


2343 


7.21 


15.79 


1.02 


0.54 


3 


-117.24 


0.45 




17 16 29.83 


+43 06 04.3 


VII-8 


7.11 


14.94 


1.10 


0.58 


4 


-117.28 


0.29 


T 


17 16 29.99 


+43 02 45.3 


1995 


8.68 


15.53 


0.96 


0.55 


4 


-124.87 


0.57 


T 


17 16 31.13 


+43 04 46.5 


3917 


7.39 


16.77 


0.93 


0.47 


2 


-123.97 


1.30 




17 16 33.36 


+43 10 21.8 


1936 


6.51 


15.48 


1.03 


0.54 


3 


-124.71 


0.54 


T 


17 16 33.89 


+43 07 04.0 


3180 


6.15 


16.37 


0.96 


0.50 


3 


-119.18 


0.72 





Note. — The full table is available on-line as a machine-readable table. 

a Number of observations. Where there is only one observation, the configuration code from Table [2] is given. For 
stars with configuration 'Z', the velocity is based on the sum of all available spectra. 



b R: Has a proper motion in iReed ( 19921 ). T: Has a proper motion in lTucholke et al.1 ( 19961 ) . v: Variable in our 
data. V: Known to be variable from another source, j: Velocities likely affected by jitter, p: Photometry is from 
a lower-quality source. 



Table 4. Data for doubtful cluster members 



11. A. 


ueci. 




r 










v 




Notes' 3 and 


(J2000.0) 


(J2000.0) 


ID 


(0 


T2 


M-T2 


C-M 


N v a 


(kms -1 ) 


(km s _1 ) 


Other Names 


(.1) 


\ z ) 


{■V 


\V 


\p) 








{a\ 
l y J 


{-\r\\ 
(10) 


I 11 ) 


17 15 59.02 


+43 09 53.3 


2171 


12.52 


15.67 


0.88 


0.55 


5 


-89.33 


0.39 


To 


17 1 R 1 Q fi/l 
1 1 ID la.041 


-\-tO 14 UO.O 


9Q1 Q 

zy iy 


10/1/1 


10. zu 


n qi 
u.y i 


u.^y 


o 
•J 


1 nn so 


l.UD 


io 


17 17 32.80 


+43 12 00.4 


3448 


6.05 


16.52 


1.07 


0.62 


3 


-94.31 


0.88 





17 18 12.18 


+42 55 04.0 


30 


17.74 


10.26 


2.03 


1.70 


1 


-108.45 


1.12 


Tpo 


17 16 29.11 


+43 09 44.9 


VI-7 


7.09 


12.66 


1.03 


0.65 


ml 


-139.57 


0.31 


RTc ZNG-4 


17 17 07.91 


+43 07 12.2 


R644 


1.00 


14.03 


1.09 


1.06 


3 


-121.84 


0.92 


Rh 


17 18 39.80 


+43 06 10.5 


1016 


17.05 


14.91 


0.70 


0.70 


J 2 


-85.97 


0.87 


Th 



^Number of observations. Where there is only one observation, the configuration code from Table is given. 



b R: Has a proper motion in iReesI ( 19921 ) . T: Has a proper motion in lTucholke et al.1 ( 19961 ). o: Outlying velocity 
c: Odd colors, h: High h\ value, p: Photometry is from a lower-quality source. 



Table 5. Data for non- members 



R.A. 


Decl. 












V 




Notes' 3 and 


(J2000.0) 


(J2000.0) 


ID 


T2 


M-T2 


C-M 




(kms -1 ) 


(kms -1 ) 


Other Names 


(!) 


(2) 


(3) 


( 4 ) 


(5) 


(6) 


(') 


(8) 


few 

(9) 


(10) 


17 14 16.43 


+43 09 40.8 


40 


10.86 


1.42 




J3a 


-65.95 


2.35 


p 


17 14 24.38 


+42 50 39.0 


594 


13.42 


1.36 




2 


-89.90 


0.44 


p 


17 14 24.44 


+42 55 45.8 


487 


13.70 


0.90 


0.47 


J3a 


-47.57 


2.27 


p 


17 14 27.31 


+42 44 50.4 


389 


13.11 


1.24 




J4 


-3.30 


1.53 


p 


17 14 30.21 


+43 03 26.5 


447 


13.43 


0.95 


0.68 


Jl 


-2.41 


0.76 


p 


17 14 32.91 


+43 01 10.6 


739 


13.90 


1.20 


0.90 


2 


-95.67 


0.96 


p 



Note. — The full table is available on-line as a machine-readable table. 



a Number of observations. Where there is only one observation, the configuration code from Table [2] is 
given. For stars with configuration 'Z', the velocity is based on the sum of all available spectra. 



m 



b R: Has a p roper motion in iReesI ( 19921 ) and probablity of membership > 90%. r : Has a proper motion 
Reesl ()1992l ) and probablity of membership < 90%. T: Has a proper motion i n iTucholke et al.l ([1996) 



and probablity of membership > 90%. t: Has a proper motion in lTucholke et al.l (|1996l ) and probablity of 
membership < 90%. v: Variable, p: Photometry is from a lower-quality source. 



c This identification is ambiguous. This star could be R137. 



Table 6. Permitted ranges for parameters 



M92 M15 



Parameter 


Models 


lower 


upper 


lower 


upper 


v (kms -1 ) 


All 


-122.58 


-119.58 


-108.5 


-105 


<7i (kms -1 ) 


All 


0.3 


12 


0.3 


20 


a 


P,PB,C,CB 


-2.5 





-2.5 





ro (') 


C,CB 


0.3 


8 


0.3 


3 


n (') 


PB,CB 


8 


14 


7.5 


15 


cr out (kms" 1 ) 


PB,CB 


0.3 


10 


0.3 


10 



Table 7. Results for Mock Data Sets 



log(p(Model|D7) /p(Best Model)/}/)) 



Model 


PL 


CPL 


PL+F 


D 


s 


P 


0.0 


-1.2 


-7.3 


-5.8 


-2.1 


C 


-1.7 


0.0 


-8.8 


-4.0 


-1.8 


PB 


-0.5 


-1.7 


-0.04 


-6.6 


-3.7 


CB 


-2.3 


-0.3 


0.0 


-5.1 


-3.4 


B N (1-20) 


-9.5 


-3.0 


-5.7 


-1.5 


0.1 


Bn (best) 


-9.9 (11) 


-3.3 (7) 


-6.1 (7) 


-1.5 (5) 


0.0 (1) 


B S (2-4) 


-10.5 


-4.3 


-6.9 


0.1 


0.0 


B s (best) 


-10.5 (4) 


-4.3 (3) 


-6.9 (4) 


0.0 (3) 


-0.1 (2) 



Table 8. Results for M92 



log 


" p{D\HI) 1 


log 


r P (h\i) i 


log 


r P {h\di) 1 




p{D\CB I) 


[p(CB\I) 


[p(CB\DI)\ 




Model 


(1) 




(2) 




(3) 



P 


-0.5 


-2.0 


-2.5 


c 


-0.2 


0.0 


-0.2 


PB 


-0.9 


-2.0 


-2.9 


CB 


0.0 


0.0 


0.0 


B 


0.1 


-1.7 


-1.6 


Bn=2 


0.0 


-2.9 


-2.9 


Bn=3 


-1.4 


-2.9 


-4.3 


Bn=4 


-1.5 


-2.9 


-4.4 


Bn=5 


-2.1 


-2.9 


-5.0 


Bn=6 


-3.2 


-2.9 


-6.1 


Bs=2 


-0.7 


-2.9 


-3.6 


Bs=3 


-1.1 


-2.9 


-4.0 


B S =A 


-1.8 


-2.9 


-4.7 



Note. — (1) Log relative likelihood, (2) Log relative prior 
probability, (3) Log relative posterior probability 
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Table 9. Mean points for M92 





Model C 






Model CB 


r 


\° V II 




(r) 




>)) 


t<r(r) 


(') 


( Irvn t; 1 ^ 


(kms 1 ) 


(kms 1 ) 


firm t; 1 ^ 


U.z4 


7.25 




1 06 




6 90 


0.93 


U.ol) 


6.54 




57 




6 46 


0.56 


1 OP. 


5.89 




^6 




5 96 


0.38 


1 no 

i.yz 


5.35 




99 

VJ. iJ 




5 48 


0.31 


O A C 


4.90 




96 

VJ.ZjVJ 




5 04 


0.28 


o.U4 


4.54 




94 

VJ. Zj'-t 




4 64 


0.26 


o.bl) 


4.24 




99 

VJ. ZjZj 




4 SI 

"-±.0 ± 


0.28 


4.1b 


3.98 




99 




S 95 


0.22 


/I TO 


3.76 




91 




O. I O 


0.23 


r: o o 

5.28 


3.58 




0.22 




3.44 


0.26 


5.84 


3.41 




0.22 




3.25 


0.24 


6.40 


3.27 




0.23 




3.07 


0.28 


6.96 


3.14 




0.24 




2.89 


0.32 


7.52 


3.03 




0.25 




2.75 


0.33 


8.08 


2.93 




0.26 




2.97 


0.72 


8.64 


2.84 




0.26 




3.32 


0.89 


9.20 


2.75 




0.27 




3.39 


0.94 


9.76 


2.68 




0.28 




3.37 


0.97 


10.32 


2.61 




0.29 




3.82 


0.98 


10.88 


2.54 




0.29 




3.81 


0.98 


11.44 


2.48 




0.30 




3.81 


0.99 


12.00 


2.43 




0.31 




3.81 


0.99 


12.56 


2.38 




0.31 




3.81 


0.99 


13.12 


2.33 




0.32 




3.82 


1.00 


13.68 


2.28 




0.32 




3.84 


1.02 


14.24 


2.24 




0.32 




3.97 


1.14 



Table 10. Results for M15 



log 


" p{D\HI) 1 


log 


r P (h\i) i 


log 


r P {h\di) 1 




p{D\CB I) 


[p(CB\I) 


[p(CB\DI)\ 




Model 


(1) 




(2) 




(3) 



P 


-0.3 


0.0 


-0.3 


c 


-0.6 


-0.7 


-1.3 


PB 


-0.3 


0.0 


-0.3 


CB 


0.0 


-0.7 


-0.7 


B 


-0.9 


-1.7 


-2.6 


Bn=2 


-3.2 


-2.9 


-6.1 


Bn=3 


-0.9 


-2.9 


-3.8 


Bn=5 


-3.0 


-2.9 


-5.9 


Bn=6 


-3.4 


-2.9 


-6.3 


Bn=7 


-2.8 


-2.9 


-5.7 


B S =2 


-1.9 


-2.9 


-4.8 


Bs=3 


-1.7 


-2.9 


-4.6 


B S =A 


-1.9 


-2.9 


-4.8 



Note. — (1) Log relative likelihood, (2) Log relative prior 
probability, (3) Log relative posterior probability 
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Table 11. Mean points for M15 





Model P 




Model PB 


r 


\° V II 




•(r) 


(a(r)) 


t<r(r) 


(') 


( Irvn t; 1 ^ 


(kms 1 ) 


( lr m t; ^ ^\ 


( lr tyi t; ^ 


n no 
U.Uo 


20.85 




4 60 


28.99 


8.38 


I). 74 


8.47 




83 


9.29 


1.08 


1.41 


6.59 




46 


6.80 


0.51 


2.07 


5.68 




39 

VJ. OA 


5.65 


0.33 




5.10 




96 

VJ.ZjXJ 


4.96 


0.26 


3.40 


4.69 




93 

VJ. S-IKJ 


4.47 


0.25 


4.06 


4.38 




0.22 


4.12 


0.25 


4.72 


4.14 




99 

U - — Z, 


3.84 


0.26 


c on 

5.39 


3.94 




93 


3.61 


0.27 


O.05 


3.77 




0.23 


3.43 


0.28 


0. f 1 


3.63 




0.23 


3.27 


0.28 


r.oo 


3.50 




0.24 


3.13 


0.29 


8.04 


3.39 




0.24 


3.05 


0.33 


8.71 


3.29 




0.25 


3.05 


0.43 


9.37 


3.20 




0.25 


3.40 


0.62 


10.03 


3.12 




0.25 


3.44 


0.65 


10.69 


3.05 




0.26 


3.65 


0.68 


11.36 


2.98 




0.26 


3.77 


0.69 


12.02 


2.92 




0.26 


3.77 


0.70 


12.68 


2.86 




0.26 


3.79 


0.70 


13.35 


2.81 




0.27 


3.79 


0.71 


14.01 


2.76 




0.27 


3.79 


0.72 


14.67 


2.71 




0.27 


3.81 


0.75 


15.34 


2.67 




0.27 


3.81 


0.75 


16.00 


2.63 




0.27 


3.81 


0.75 


16.66 


2.59 




0.27 


3.81 


0.75 



